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Abstract 

This  research  examined  the  theory  and  application  of  using  orthogonal  fre¬ 
quency  division  multiplexing  (OFDM),  or  discrete  multi-tone  (DMT),  frequency- 
domain  equalization  (FEQ)  with  two  communications  systems  that  inherently  possess 
unused  carrier  frequencies,  or  null-tones,  in  their  respective  transmission  frequencies. 
The  fundamental  DMT-FEQ  property  relies  on  null-tones  to  equalize  a  non-ideal 
channel  and  mitigate  the  effects  of  interchannel  interference  (ICI),  intersymbol  in¬ 
terference  (ISI),  and  noise.  The  two  communications  systems  investigated  were  a 
Hadamard  encoded  code  division  multiple  access  (CDMA)  communications  system 
with  up  to  32  synchronous  users  and  a  transform  domain  communications  system 
(TDCS)  with  only  one  user.  Both  communications  systems  were  simulated  while 
operating  with  real  channel  data  corrupted  by  noise. 

Simulation  results  showed  that  the  Hadamard  encoded  CDMA  system  worked 
well  with  DMT-FEQ  only  when  the  Hadamard  code  set  was  used  to  construct  a 
transmission  signal  that  obeyed  DMT-FEQ  null-tone  conditions  in  conjunction  with 
a  vector  estimation  method.  Simulation  results  also  showed  that  a  TDCS  using  tradi¬ 
tional  pseudo-random  phase  component,  and  traditional  spectral  mask  with  consecu¬ 
tive  null-tones,  did  not  work  well  with  DMT-FEQ.  Modifications  to  the  TDCS  model 
revealed  that  a  TDCS  with  a  conjugate-symmetric  phase  component  in  conjunction 
with  a  modified  spectral  mask  with  consecutive  null-tones  and  forced  null-tones  pro¬ 
vided  acceptable  results  when  equalizing  with  DMT-FEQ.  The  DMT-FEQ  may  be 
suitable  for  covert  applications,  such  as  TDCS,  when  modifications  to  TDCS’  phase 
component  and  forced  null-tones  in  its  spectral  mask  are  made. 
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Applying  Frequency-Domain  Equalization 


TO 

Code-Division  Multiple  Access 
and 

Transform-Domain  Communications  Systems 

I.  Introduction 

1.1  Background 

The  discrete-multitone  (DMT)  frequency-domain  equalizer  (FEQ)  was  intro¬ 
duced  as  an  alternative  to  existing  asymmetric  digital  subscriber  line  (ADSL)  equal¬ 
ization  techniques  in  order  to  decrease  the  latency  effects  of  equalization  in  the  time- 
domain  [20-22],  The  DMT-FEQ  is  not  classified  as  a  blind  equalizer  since  it  relies  on  a 
channel  identihcation  and  used  carrier  frequencies  and  null-tones.  The  DMT-FEQ  was 
implemented  in  this  research  effort  due  to  its  authors’  claims  of  its  ability  to  perfectly 
equalize  a  channel  to  remove  interchannel  interference  (ICI)  and  intersymbol  inter¬ 
ference  (ISI)  in  a  communications  system  that  possesses  unused  carrier  frequencies, 
or  null-tones.  Two  communications  systems  that  inherently  possess  null-tones  are  a 
Hadamard  encoded  code-division  multiple  access  (CDMA)  communications  systems 
and  transform- domain  communications  systems  (TDCS). 

CDMA  techniques  were  first  introduced  in  military  applications  during  World 
War  II  [14],  Currently,  CDMA  is  used  by  both  military  and  civilian  applications 
around  the  globe.  TDCS  is  designed  for  military  applications  that  require  low  prob¬ 
ability  of  detection  (LPD)  and  low  probability  of  intercept  (LPI). 
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1.2  Problem  Statement 


The  primary  objective  of  this  research  effort  was  to  combine  a  Hadamard  en¬ 
coded  CDMA  communications  system  with  the  DMT-FEQ  described  in  [20-22].  The 
secondary  objective  was  the  use  of  DMT-FEQ  in  LPD/LPI  military  applications,  such 
as  TDCS. 

1.3  Assumptions 

Three  major  assumptions  were  made  in  order  to  progress  through  this  research 
in  a  realistic  and  timely  fashion:  the  channel  data  that  was  used  was  time  invariant 
and  was  also  shortened  to  maintain  compatibility  with  the  DMT-FEQ  (described  in 
section  3.1.1),  the  noise  inserted  into  the  receiver  in  the  simulations  was  Gaussian,  and 
the  noise  used  in  the  simulations  involving  the  first  half  of  the  Hadamard  spreading 
code  set,  hi:32,  was  assumed  to  be  zero-mean  with  unit  variance  in  order  to  construct 
an  invertible  noise  covariance  matrix.  Also,  the  channel  was  assumed  to  be  known 
exactly. 

1.4  Methodology 

The  methodology  employed  in  this  research  effort  explored  two  communications 
systems,  Hadamard  encoded  CDMA  and  TDCS,  that  inherently  possess  null-tones  in 
their  respective  transmission  carrier  frequencies,  used  in  conjunction  with  the  DMT- 
FEQ.  With  respect  to  the  Hadamard  encoded  CDMA  system,  intricacies  of  the  trans¬ 
mitter,  properties  of  the  DMT-FEQ,  impact  of  equalizing  received  signals  with  the 
DMT-FEQ,  derivations  of  both  the  scalar  and  vector  maximum  likelihood  estimates 
(MLE),  and  the  requirement  for  a  post-estimation  decision  function  were  all  thor¬ 
oughly  explored.  And,  with  respect  to  the  TDCS  model,  incremental  changes  in  the 
spectral  mask  and  the  phase  component  was  conducted  in  order  to  try  to  identify 
non-compatibility  issues  with  the  DMT-FEQ  and  a  traditional  TDCS  system. 
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1 . 5  Materials 


The  following  materials  were  vital  to  this  research  effort:  Dell  Precision  380 
desktop  computer  with  MATLAB  version  7.4.0.287  (R2007a)  for  simulations. 

1 . 6  Overview 

Chapter  one  is  the  introduction.  Chapter  two  reviews  existing,  relevant  ma¬ 
terial  to  the  research  presented  in  this  thesis  and  is  intended  to  provide  the  reader 
with  a  thorough  frame  of  reference  for  the  proceeding  material.  Chapter  three  de¬ 
scribes  the  implementation  of  the  DMT-FEQ  in  conjunction  with  two  different  types 
of  communications  system  models  that  inherently  possess  null-tones  in  their  respec¬ 
tive  transmission  carrier  frequencies.  Chapter  four  describes  simulation  results  of 
equalizing  a  Hadamard  encoded  CDMA  communications  system  and  a  TDCS  with 
the  DMT-FEQ.  Chapter  five  concludes  this  research  effort. 
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II.  Literature  Review 


The  research  presented  in  this  thesis  involves  the  discrete-mnltitone  (DMT), 
also  known  as  orthogonal  frequency  division  multiplexing  (OFDM)  frequency- 
domain  equalizer  (FEQ)  used  in  conjunction  with  communication  systems  that  in¬ 
herently  possess  null-tones  in  their  respective  transmission  signals.  Therefore,  in 
order  to  provide  a  solid  frame  of  reference  for  the  research  presented  in  this  the¬ 
sis,  the  following  areas  of  signal  processing  and  communications  will  be  briefly  re¬ 
viewed:  DMT-FEQ,  code-division  multiple  access  (CDMA)  communications  systems, 
quadrature  phase  shift  keying  (QPSK),  quadrature  amplitude  modulation  (QAM), 
and  transform- domain  communications  systems  (TDCS). 

2.1  Discrete- Multitone  Frequency- Domain  Equalization 

The  communications  systems  used  in  the  simulations  were  modeled  with  the 
DMT-FEQ  as  described  in  [20-22],  The  DMT-FEQ  achieves  signal  reconstruction 
by  utilizing  samples  from  normally  unused  carriers,  or  null-tones,  to  equalize  inter¬ 
channel  interference  (ICI)  and  intersymbol  interference  (ISI).  The  DMT-FEQ  is  not 
classified  as  a  blind  equalizer  since  it  relies  on  a  channel  identification  and  used  carrier 
frequencies  and  null-tones. 

DMT-FEQ  signal  reconstruction  can  be  achieved  as  long  as  either  of  the  two- 
case  null-tone  rules  and  the  equidistant  null-tone  rule  noted  in  [20-22]  are  not  violated. 
These  null-tone  rules  state  that  case  1  is  K  >  Lc  —  1,  Lc  —  1  <  2N,  case  2  is  K  > 
2 N,  Lc  —  1  >  2N,  and  the  equidistant  null-tone  rule  requires  —  1  used  carrier 
frequencies  between  each  null-tone,  where  K  denotes  the  total  number  of  null-tones 
present,  Lc  is  the  length  of  the  channel  impulse  response  (CIR),  and  N  is  the  total 
number  of  the  communications  system’s  used  carrier  frequencies. 

The  simulation  results  noted  in  [20-22]  primarily  focused  on  an  OFDM  system’s 
performance  with  respect  to  its  data  rate  (Kbps)  versus  its  introduced  redundancy 
(^y),  where  CP  denotes  the  length  of  the  cyclic  prefix.  The  simulations  noted  in 
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Figure  2.1:  TDMA  users  are  restricted  to  individual  time  slots 
to  transmit  and  receive  data.  Typically,  guard  time  (GT)  is  used 
in  TDMA  systems  to  prevent  inter-user  interference.  FDMA 
users  are  restricted  to  individual  carrier  frequencies,  fc,  to  trans¬ 
mit  and  receive  data. 


[20-22]  maximized  an  OFDM  system’s  data  rates  for  a  given  probability  of  bit  errors, 
Pb.  See  Fig.  8  of  [21]. 


2.2  Code-Division  Multiple  Access  Communications  Systems 

Communications  systems  that  employ  CDMA  are  also  referred  to  as  spread  spec¬ 
trum  multiple  access  communications  systems  due  to  their  ability  of  allowing  multiple 
signals  to  simultaneously  occupy  the  same  radio  frequency  (RF)  bandwidth  without 
interfering  with  one  another  [14,19].  CDMA  does  not  have  restrictions  on  time  or 
bandwidth  that  time-division  multiple  access  (TDMA)  or  frequency-division  multiple 
access  (FDMA),  respectively  possess.  As  shown  in  Fig.  2.1,  users  of  TDMA  com¬ 
munications  systems  and  FDMA  communications  systems  are  restricted  on  time  and 
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Figure  2.2:  CDMA  signal  detection,  (a)  Wideband  spectrum 
into  the  receiver  consists  of  wanted  and  unwanted  signals  that 
are  each  spread  by  their  unique  code  with  code  rate,  Rch-  (b) 
Narrowband  spectrum  in  the  receiver  after  correlation  with  the 
user’s  unique  code  occupies  a  bandwidth  centered  at  an  inter¬ 
mediate  frequency  (IF),  which  is  then  applied  to  a  demodulator. 


frequency  allocation,  respectively.  CDMA’s  RF  spectrum  at  the  input  to  the  receiver 
and  the  RF  spectrum  after  correlation  with  the  correct  synchronized  pseudorandom 
number  (PN)  code  is  shown  in  Fig.  2.2. 

2.2.1  Direct  Sequence  Code  Division  Multiple  Access.  A  typical  direct 
sequence-CDMA  (DS-CDMA)  system  is  shown  in  Fig.  2.3,  where  J  is  the  total  num¬ 
ber  of  users.  Fig.  2.3  shows  user  l’s  data  modulated,  ui(t),  then  multiplied  with  a 
unique  spreading  signal,  hi(t).  Users  2  through  J  may  also  simultaneously  spread 
their  modulated  data  with  spreading  codes  unique  to  them.  As  shown  in  Fig.  2.2, 
multiplication  in  the  time-domain,  u(t)h(t),  results  in  convolution  in  the  frequency- 
domain  [18],  therefore,  assuming  u(t)  relatively  occupies  a  narrow  portion  of  the 
frequency  spectrum  compared  with  h(t),  the  product  signal,  u(t)h(t),  will  be  approx¬ 
imately  the  bandwidth  of  h(t). 

In  the  absence  of  noise,  the  signal  present  at  the  receiver,  y(t),  is  defined  as 

j 

y(t)  =  2Jun(t)hn(t).  (2.1) 

n=l 
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\Thj(t)  hj(t)  dt= 

Jo  '  1 0  for  it) 


Carrier: 

A  cos  u>ot 

I 

Modulator 

nr 

Transmit 

Data 


uj(t)  hj(t) 


Modulated  signal: 

Ui(t)  =  Ai(t)  cos  [uj0t  +  <Pi(t)] 


<*> 

t 


u2(t)  h2(t) 
ui(t)  h 


I  .  • 

H> 


i — -<S> 


ui(t)  hi2(t) 
u2(t)  h2(t)  h-i(t) 

Uj(t)  hj(t)  hrf) 


To  Demodulator 


Spreading  Code:  Despreading  Code: 

hi(t)  htf) 


Transmitter 

Receiver 

Figure  2.3:  A  direct  sequence  code  division  multiple  access 

(DS-CDMA)  system.  Shown  for  spreading/despreading  user  l’s 
transmission. 


The  first  stage  of  the  receiver  in  Fig.  2.3  multiplies  the  incoming  signal,  (2.1), 
with  a  unique  despreading  code,  h\(t),  for  example,  to  obtain  the  desired  signal, 


plus  unwanted  signals, 


(2.2) 


u2(t)h2{t)h1(t)  +  u3(t)h3(t)hi(t)  H - b  uj(t)hj(t)hi(t).  (2.3) 

If  the  spreading/despreading  codes  are  chosen  with  orthogonal  properties,  as 
described  in  section  2.2.2,  the  desired  signals  can  be  extracted  perfectly  in  the  absence 
of  noise.  This  is  because  JQT  hf(t)dt  =  1  and  J0T  hi{t)hj{t)dt  =  0  for  i  ^  j  [14], 

The  frequency-domain  view  of  the  DS-CDMA  receiver  is  shown  in  Fig.  2.2. 
Fig.  2.2(a)  depicts  the  wideband  input  into  the  receiver  that  consists  of  wanted  and  un¬ 
wanted  signals  that  are  each  spread  by  their  own  code  with  code  rate,  Rch ■  Fig.  2.2(b) 
illustrates  the  RF  spectrum  after  correlation  with  the  user’s  unique  code.  The  de¬ 
sired  signal  in  (2.2)  is  centered  at  an  intermediate  frequency  (IF)  and  occupies  the 
information  bandwidth.  Typically  the  signal  produced  by  (2.2)  is  then  applied  to 
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a  demodulator  with  its  bandwidth  equal  to  the  information  bandwidth  to  pass  the 
desired  despread  signal.  The  undesired  signals  of  (2.3)  remain  spread  by  hi(t)hn(t) 
and  only  the  portion  of  the  spectrum  of  the  unwanted  signals  falling  in  the  informa¬ 
tion  bandwidth  of  the  receiver  will  cause  interference  with  the  desired  signal  in  the 
demodulation  process  [14]. 

2.2.2  Orthogonal  Signals.  An  example  of  orthogonal  signals  is  shown  in 
Fig.  2.4,  where  u\(t)  and  u2{t)  are  constructed  from  pulse  waveforms  described  by 
ui(t)  —  p{t)i  0  <  t  <  T  and  u2{t)  =  pit  —  |-),0<t<T,  Tis  the  symbol  period,  and 
p(t)  is  a  pulse  with  duration  r  =  As  described  in  [14,18],  a  set  of  equal  energy 
signals  Ui(t),  where  *  =  1,  2, . . . ,  J,  is  orthogonal  if  and  only  if 


(2.4) 


where  is  the  cross-correlation  coefficient,  and  E  is  the  signal  energy,  defined  as 


(2.5) 


Jo 

An  example  of  computing  the  cross-correlation  between  two  orthogonal  signals 
at  the  waveform  level  is  illustrated  in  Fig.  2.5(a). 

2.2.3  Waveform  Coding.  The  ultimate  goal  of  waveform  coding  is  to  produce 
the  smallest  cross-correlation  result  from  (2.4)  among  all  pairs  of  signals.  According 
to  [14,18,23],  it  is  possible  to  make  the  cross-correlation  coefficient  equal  to  zero. 

(2.4)  presents  a  definition  of  orthogonality  in  terms  of  waveforms  iq(t)  and  Uj(t), 
where  i,  j  =  1, . . . ,  M,  and  M  is  the  size  of  the  waveform  set.  Each  waveform  in  the 
set  {ui(t)}  consists  of  a  sequence  of  pulses,  where  each  pulse  is  designated  with  a  level 
+  1  or  -1,  which  in  turn  represents  the  binary  digit  1  or  0,  respectively.  When  the  set 
is  expressed  in  binary,  (2.4)  becomes 


Analytic  Waveform  Vector 

Representation  Representation  Representation 


ui(t)  =  P(t) 


u2(t)  =  p(t  -  L) 

0<t<T  2 

(a) 

Figure  2.4:  Binary  Orthogonal  Signals,  (a)  Analytic  repre¬ 

sentation  of  orthogonal  signals,  Ui(t)  and  w2(t).  (b)  Ui(t)  and 
u2(t)  cannot  interfere  with  one  another  because  they  are  disjoint 
in  time,  (c)  Vector  representation  illustrating  the  perpendicular 
relationship  between  signals  ui(t)  and  112(f) . 


ui(t) 


time 


T  T 
2 

u2(t) 


time 


T  T 
2 

(b) 


(c) 


jj  digit  agreements  —  |j  digit  disagreements 
13  fi  total  digits  in  sequence 

An  example  of  computing  the  cross-correlation  between  two  orthogonal  signals 
at  the  bit  level  is  illustrated  in  Fig.  2.5(b). 

2.2.4  Orthogonal  Codes.  As  described  in  [14],  orthogonal  codewords  of  two 
digits  each  can  be  used  to  transform  a  one-bit  data  set,  Bi,  described  by  the  rows  of 
the  matrix  Hx  as 


1,  i  =  j 
0.  » 7 ij 


(2.6) 


1 

O 

_ 1 

1 

O 

1 

O 

Bi  = 

1 

,  H1  = 

0 

1 

(2.7) 


(2.6)  can  be  used  to  verify  the  orthogonality  of  the  codeword  set,  Hi.  To  encode  a 
2-bit  data  set,  B2.  H2  will  be  defined  as 
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“Waveform”  Level  Calculation 


‘Bit”  Level  Calculation 


Vi  =  [0  0  1  1] 

v2  =  [0  1  1  0] 
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V  Bi2Ts  ■  B2Ts  ; 
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(a)  (b) 

Figure  2.5:  Cross-Correlation  Between  Orthogonal  Signals, 

(a)  Waveform  level  calculation  of  the  cross-correlation  between 
two  orthogonal  signals,  (b)  Bit  level  calculation  of  the  cross¬ 
correlation  between  two  orthogonal  signals. 
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H2’s  lower  right  quadrant  is  the  complement  of  the  prior  codeword  set.  Con¬ 
tinuing  the  construction  rule  to  obtain  an  orthogonal  set  H3  for  a  3-bit  data  set,  B3, 
produces 
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(2.9) 


In  general,  construction  of  a  codeword  set,  H„,  of  dimension  2n  x  2n,  is  called 
a  Hadamard  matrix.  For  a  n-bit  data  set,  Hn  is  defined  as 


H 


n 


H 


n— 1 


H 


71  —  1 


(2.10) 


Hn_i 

Each  pair  of  words  in  each  codeword  set  Hi,  H2,  . . .,  Hn  has  equal  number  of  digit 
agreements  and  digit  disagreements  [7],  and,  as  defined  in  (2.6),  ctl  =  0  for  i  ^  j, 
which  makes  each  of  the  codeword  sets  orthogonal. 


2.3  Quadrature  Phase  Shift  Keying 

Consider  an  original  data  stream,  d0l  d\,d2, . . .  consisting  of  bipolar  pulses  shown 
in  Fig.  2.6(a);  that  is,  the  values  of  d(t)  are  +1  and  —1,  representing  binary  one  and 
zero,  respectively.  d{t)  can  be  divided  into  an  in-phase  stream,  d/(f),  and  a  quadrature 
stream,  dg(t),  illustrated  in  Fig.  2.6(b),  as 
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Figure  2.6:  Illustrations  of  BPSK.  (a)  BPSK  original  data 

stream  consisting  of  bipolar  pulses,  (b)  Pulse  stream  divided 
into  an  in-phase  stream,  di(t),  and  a  quadrature  stream,  d,Q(t). 
(c)  A  BPSK  modulator. 


di(t)  =  d0,d2,d4,  •  •  • 


(2.11) 


dg(t)  —  di,  d3,  d5, . . . , 


where  dj(t)  is  the  even  bits  and  dQ[t)  is  the  odd  bits  of  the  original  data  stream,  d{t). 

A  convenient  orthogonal  realization  of  a  QPSK  waveform,  u(t),  is  achieved  by 
amplitude  modulating  the  in-phase  and  quadrature  data  streams  onto  the  cosine  and 
sine  functions  of  a  carrier  wave,  as  defined  in  [14]  as 

u(t)  =  ~^=  dj(t )  cos  (u0t  +  0  +  -^=  dQ(t)  sin  +  0  .  (2-12) 
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TT  , 


sin  (u)0t  +-j) 


Figure  2.7:  QPSK  and  OQPSK  Signal  Space.  Due  to  the 

fact  that  cos  (a }0t  +  jj  and  sin  (a j0t  +  |)  are  orthogonal,  the  two 
BPSK  signals  can  be  detected  separately. 


Using  trigonometric  identities  in  [13,15],  (2.12)  can  be  written  as 


u(t)  =  cos  [u0t  +  0(f)] ,  (2-13) 

which  is  shown  in  in  Fig.  2.6(c). 

The  pulse  stream,  di(t),  amplitude- modulates  the  cosine  function  with  an  am¬ 
plitude  of  ±1,  which  is  equivalent  to  a  binary  phase  shift  keying  (BPSK)  waveform 
because  the  phase  of  the  cosine  function  is  shifted  by  either  0  or  n.  Similarly,  the 
pulse  stream,  dqit),  modulates  the  sine  function  to  produce  a  BPSK  waveform  that 
is  orthogonal  to  the  cosine  function.  The  summation  of  these  two  orthogonal  compo¬ 
nents  of  the  carrier  yeilds  the  QPSK  waveform.  The  value  of  9(t)  will  correspond  to 
one  of  the  four  possible  combinations  of  dj(t)  and  dq{t)  in  (2.12):  6{t )  =  0,  ±|,  or  7r. 
The  signal  vectors  are  shown  in  Fig  2.7. 
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(a)  (b) 

Figure  2.8:  Illustrations  of  QAM.  (a)  4-ary  signal  space.  Of¬ 
tentimes  for  convenience,  4-QAM  data  is  represented  by  ±l±j. 

(b)  A  general  QAM  modulator. 

2-4  Quadrature  Amplitude  Modulation 

QAM  consists  of  two  independently  amplitude-modulated  carriers  in  quadra¬ 
ture,  which  is  considered  a  logical  extension  of  QPSK  [14].  Each  block  of  m  bits  (m 
assumed  even)  can  be  split,  into  two  (y)-bit  blocks  which  use  (y)-bit  digital-to-analog 
converters  to  provide  the  required  modulating  voltages  for  the  carriers.  At  the  re¬ 
ceiver,  each  of  the  two  signals  is  independently  detected  using  matched  filters  [14]. 
Fig.  2.8(a)  illustrates  a  two-dimensional  signal  space  and  a  set  of  4-ary  QAM  signal 
vectors  or  points  arranged  in  a  rectangular  constellation.  Oftentimes,  4-QAM  data 
is  conveniently  represented  as  ±1  ±  j.  A  general  QAM  modulator  is  shown  in  Fig. 
2.8(b). 


2.5  Transform  Domain  Communication  System 

TDCS  is  a  communications  system  that  is  designed  for  covert  operations  be¬ 
cause  of  the  security  provided  to  its  users  by  its  inherent  low  probability  of  detection 
(LPD)  and  low  probability  of  intercept  (LPI)  [2,3].  TDCS  is  a  LPI/LPD  communica¬ 
tions  system  because  of  its  ability  to  utilize  a  noise-like  time-domain  communications 
waveform  while  simultaneously  avoiding  interfering  signals  in  the  RF  spectrum  of 
interest . 
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Figure  2.9:  Typical  configuration  of  a  TDCS  transmitter. 


2.5.1  TDCS  Transmitter.  A  typical  TDCS  transmitter  is  illustrated  in  Fig. 
2.9  [3],  where  xk  represents  a  TDCS  user’s  binary  data  to  transmit,  uk  denotes  a 
TDCS  user’s  modulated  data  symbol,  and  v(t)  represents  the  TDCS  transmission 
signal. 


2.5. 1.1  Spectrum  Estimate.  There  are  several  methods  of  obtaining 
information  on  the  frequency  spectrum  available  for  the  TDCS  transmissions  [2,6,12, 
16],  but  basically  the  regions  of  the  frequency  spectrum  containing  interference  are 
identified  to  be  avoided. 

2.5. 1.2  Spectrum  Magnitude.  The  purpose  of  the  spectrum  magni¬ 
tude  is  to  construct  a  frequency-domain  communications  waveform,  referred  to  by 
Chakravarthy  et  al.  as  the  fundamental  modulation  waveform  (FMW)  [2].  This 
identifies  regions  of  the  frequency  spectrum  to  avoid  based  on  a  the  spectrum  esti¬ 
mate,  which  yields  a  communications  waveform  whose  power  spectral  density  (PSD) 
contains  zero  power  in  the  portion  of  the  spectrum  identified  as  having  interference 
present.  This  produces  a  vector  of  magnitude  components,  A'{uS),  valued  at  either 
1  or  0,  whose  values  correspond  to  interference  at  given  frequencies  c a  based  on  the 
spectrum  estimate  that  are  below  or  above,  respectively,  a  specific  LPD/LPI  thresh- 
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Figure  2.10:  An  example  of  TDCS’  functionality  of  the  spec¬ 
trum  estimate  and  spectrum  magnitude. 

old  [3].  The  spectrum  estimate  and  spectrum  magnitude  are  illustrated  in  Fig.  2.10. 
Note  that  this  creates  a  FMW  with  the  null-tones  required  by  the  DMT-FEQ. 


2.5. 1.3  Pseudo-Random  Phase.  This  provides  the  covert  applicability 
to  TDCS  via  a  complex  pseudorandom  phase  vector  because  this  produces  a  noise¬ 
like  time-domain  communications  waveform  [2].  As  described  in  [9],  TDCS  multiple 
access  networks  can  then  be  constructed  by  using  the  pseudorandom  phase  vector  as  a 
spreading  code.  The  pseudorandom  phase  vector  must  be  generated  on  a  symbol-by¬ 
symbol  basis  in  order  to  ensure  TDCS  transmissions  maintain  a  noise-like  appearance 
over  time,  which  secures  TDCS’  LPI/LPD  status. 

2. 5. 1.4  Scaling.  Scaling  ensures  an  appropriate  amount  of  energy  is 
contained  in  the  communications  symbol  to  be  transmitted,  and  that  all  communica¬ 
tions  symbols  contain  uniform  energy  [2].  The  uniform  energy  throughout  a  TDCS 
waveform  is  important  because  of  the  noise-like  structure  of  a  TDCS  waveform. 
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Figure  2.11:  Typical  configuration  of  a  TDCS  receiver.  Note 
that  the  devices  enclosed  by  the  dashed  line  perform  FMW  gen¬ 
eration  identical  to  the  process  used  by  the  TDCS  transmitter. 


2.5. 1.5  IDFT  and  Modulator.  The  purpose  of  the  inverse  discrete 
Fourier  transform  (IDFT)  and  modulator  are  self-explanitory,  and  it  is  assumed  the 
reader  need  no  further  explanations  on  these  functions. 

2.5.2  TDCS  Receiver.  A  typical  TDCS  receiver  is  illustrated  in  Fig.  2.11  [3], 
where  Uk  denotes  the  estimate  of  a  TDCS  user’s  modulated  data  symbol  and  Xk 
represents  the  estimate  of  a  TDCS  user’s  binary  data  that  was  received. 
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III.  Discrete-Multitone  Frequency-Domain  Equalization 


Implementation 

his  chapter  describes  the  implementation  of  the  discrete-multitone  (DMT) 


_L.  frequency-domain  equalizer  (FEQ)  described  in  [20-22]  in  conjunction  with  two 
different  types  of  communications  systems  that  inherently  possess  null-tones  in  their 
respective  transmission  carrier  frequencies.  The  first  system  is  a  code-division  multiple 
access  (CDMA)  system  utilizing  Hadamard  spreading  codes  for  one  to  32  synchronous 
users  and  the  second  system  is  a  transform-domain  communications  system  (TDCS) 
for  one  user. 

The  Hadamard  encoded  system  is  illustrated  in  Fig.  3.1  and  consists  of  system 
users’  binary  data,  4-QAM  modulator,  channel,  DMT-FEQ,  post-equalized  null-tone 
insertion,  maximum  likelihood  estimation  (MLE),  decision  function,  and  a  4-QAM 
demodulator.  Note  that  the  MLE  has  three  possible  estimators:  1)  a  scalar  MLE 
(described  in  section  3.1.4)  which  requires,  2)  an  initial  scalar  approximate  MLE 
(described  in  section  3. 1.4. 4),  or  3)  a  vector  MLE  (described  in  section  3.1.5). 

The  TDCS  implementation  is  illustrated  in  Fig.  3.14  and  consists  of  a  transmit¬ 
ter  (described  in  section  2.5.1),  DMT-FEQ,  and  a  receiver  (described  in  section  2.5.2). 

3.1  Hadamard  Encoded  CDMA  System  with  DMT-FEQ 

The  DMT-FEQ  described  in  [20-22]  completely  eliminates  interchannel  inter¬ 
ference  (ICI)  and  intersymbol  interference  (ISI)  when  N  carrier  frequencies  are  used 
and  K  unused  carrier  frequencies  (or  null  tones)  are  present.  The  total  number  of 
carrier  frequencies  is  M  and  is  defined  as  M  =  N  +  K.  Note  that  certain  conditions 
must  be  met  in  order  for  the  DMT-FEQ  to  mitigate  ICI  and  ISI,  which  are  discussed 
in  section  3.1.1. 

The  spreading  codes  used  to  distinguish  individual  users  are  assumed  to  be 
synchronous  and  can  therefore  be  modeled  as  an  orthogonal  Hadamard  matrix  [8]. 
As  shown  in  Fig.  3.2,  Hadamard  codes  inherently  possess  null-tones  in  the  frequency 
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Figure  3.1:  A  complete  view  of  the  Hadamard  encoded  CDMA 
communications  system.  Note  that  when  the  post-equalized  re¬ 
ceived  signal,  y,  enters  the  estimator,  it  is  either  routed  through 
the  scalar  estimator  or  the  vector  estimator.  Whichever  estima¬ 
tor  the  equalized  signal  is  routed  through,  the  output  is  u ml- 


domain  which  make  them  an  ideal  candidate  as  a  method  to:  1)  distinguish  data 
between  users  in  a  multi-user  system  and  2)  keep  null-tones  intact  in  order  to  ex¬ 
ploit  properties  of  the  DMT-FEQ  to  mitigate  1S1  and  ICI  from  a  non-ideal  channel 
corrupted  by  noise.  More  discussion  of  correctly  chosen  combinations  of  Hadamard 
codes  will  be  available  in  section  3.1.1. 
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Hadamard  Hadamard 

code  1  code  M 


I  I  I  I  I  I  I  I  I  I  I  I 
I  I  I  I  I  I  I  I  I  I  K> 


1.  i  I.  1  1  1.  1. 


I  I  I  I  I  I  I  I  1 


d3  =  nonzero  data 

Figure  3.2:  A  view  of  Hadamard  codes  in  the  frequency  do¬ 
main  show  their  inherent  null  tones.  Shown  for  M  =  16. 

3.1.1  Transmitter.  The  time-domain  transmission  signal,  denoted  by  the 
vector,  v,  is  the  Hadamard  encoded  multi-user  4-QAM  vector  prior  to  transmission 
and  is  defined  as 


j 

v  =  hn,  (3.1) 

n=  1 

where  J  is  the  total  number  of  users  in  the  system,  un  is  the  n-th  user’s  4-QAM 
symbol,  and  hn  is  n-th  user’s  Hadamard  codeword.  (3.1)  is  similar  to  (2.1)  in  all 
respects,  except  (3.1)  uses  vector  notation. 

It  is  important  to  correctly  create  v  so  that  the  DMT-FEQ  will  be  able  to 
equalize  the  received  signal  properly.  A  correctly  chosen  combination  of  Hadamard 
codes  is  one  where  equation  (3.1)  keeps  an  appropriate  number  of  null-tones  intact 
in  the  frequency-domain  of  the  transmission  vector  that  does  not  violate  either  of 
the  two-case  null-tone  rules  noted  in  [20-22]:  Case  one  is  K  >  Lc  —  1,  Lc  —  1  <  2IV, 
Case  two  is  K  >  2 N,  Lc  —  1  >  2IV,  and  Lc  is  the  length  of  the  channel  impulse 
response  (CIR).  In  other  words,  case  one  states  the  number  of  null-tones  must  be 
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Last  Half  of  Hadamard 
Codes  Used 


J=1,  2,  4,  8,  16,  32 
N=32 
K=32 


m  =  Used  Carrier  Frequency 

(a)  (b) 


Figure  3.3:  A  frequency-domain  view  of  users’  null-tone  im¬ 
pact  in  a  Hadamard  encoded  CDMA  communications  system. 
Shown  for  M  =  64.  (a)  Null-tones  vary  when  spreading  codes 
from  the  first  half  of  a  Hadamard  matrix  are  used,  (b)  Null-tones 
remain  uniform  as  long  as  the  users,  J  <  32. 


equal  to,  or  greater  than,  the  length  of  the  CIR  minus  1,  and  the  length  of  the  CIR 
minus  one  must  be  less  than  twice  the  number  of  used  carrier  frequencies.  Case  two 
states  the  number  of  null-tones  must  greater  than,  or  equal  to,  twice  the  amount  of 
used  carrier  frequencies  and  the  CIR  minus  one  must  be  greater  than,  or  equal  to, 
twice  the  amount  of  used  carrier  frequencies.  Case  two  only  applies  when  N  <  4^-, 
which  is  an  impractical  method  of  utilizing  the  allocated  frequency  spectrum  of  a 
communications  system  if  only  |  of  the  total  carrier  frequencies  are  utilized.  Case 
two  will  only  apply  when  spreading  codes  from  the  first  half  of  the  Hadamard  matrix 
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J{v} 


j{v}=Eu„-m} 


Ug  © 


+ 


+ 


K=  8 

=  nonzero  data 


User  8 


Figure  3.4:  A  frequency-domain  view  of  a  transmission  vector, 
JF{v},  created  by  a  correctly  chosen  combination  of  Hadamard 
codes  for  a  CDMA  communication  system  with  J  —  8  users, 
Lc  =  9,  M  =  16  carrier  frequencies,  N  =  8  used  carrier  fre¬ 
quencies,  and  K  =  8  null-tones.  This  is  a  correct  combination 
of  Hadamard  codes  because  it  does  not  violate  DMT-FEQ  re¬ 
quirements  K  >  Lc  —  1  and  Lc  —  1  <  2N . 


are  used  and  the  users,  J  <  16.  Please  refer  to  Fig.  3.3  for  an  illustration  on  the 
impact  of  the  number  of  system  users  on  the  frequency-domain  characteristics  of  a 
Hadamard  encoded  CDMA  communications  system’s  transmission  signals. 

Proceeding  with  a  Hadamard  encoded  CDMA  communications  system  that  fol¬ 
lows  the  case  one  rule  set,  two  options  exist  for  a  transmission  vector,  v,  composed  of 
correctly  chosen  combinations  of  Hadamard  codes  shown  in  Fig.  3.2,  where  Lc  =  9  and 
J  =  8  users.  Note  that  the  numbers  provided  are  a  simplified  example  for  illustration 
purposes  only.  Since  M  =  16,  one  option  is  to  use  the  first  half  of  the  Hadamard  ma¬ 
trix  (hi:8)  which  will  make  (3.1)  become  v  =  Yln=i  urMn  =  «ihi  +  U2I12  +  •  — V  Mshs, 
and  the  second  option  is  to  use  h9:i6  which  will  make  (3.1)  become  v  =  Yhn=\  unhn  = 
-Uihg  +  -^hio  +  •  — b  rtshie-  Both  of  these  options  do  not  violate  case  1  because  both 
instances  of  v  contain  eight  null-tones,  that  is  K  =  8  in  both  examples  given.  Fig.  3.4 
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J  I 

J{v}=Zu„-J{hn}  df> 


J{v} 


K=1 

=  nonzero  data 


Figure  3.5:  A  frequency-domain  view  of  a  transmission  vec¬ 

tor,  JF{v},  created  by  an  incorrectly  chosen  combination  of 
Hadamard  codes  for  a  CDMA  communication  system  with  J  =  8 
users,  Lc  =  9,  M  =  16  carrier  frequencies,  N  =  9  used  car¬ 
rier  frequencies,  and  K  =  7  null-tones.  This  is  an  incorrect 
combination  of  Hadamard  codes  because  it  violates  DMT-FEQ 
requirements  K  >  Lc  —  1  and  Lc  —  1  <  2N. 

shows  v  constructed  when  Hadamard  codes  nine  through  sixteen  are  used  from  a 
16  x  16  Hadamard  matrix. 

From  the  actual  64  x  64  Hadamard  matrix  used  in  the  simulated  system  model, 
Fig.  3.3(a)  and  Fig.  3.3(b)  depicts  the  following  two  scenarios  that  apply  to  the  case 
one  null-tone  rule:  1)  J  =  32  users  when  the  first  half  of  the  Hadamard  codes  are 
used  and  2)  J  <  32  when  the  last  half  of  the  Hadamard  codes  are  used.  Fig.  3.3(a) 
illustrates  when  the  case  two  null-tone  rule  applies:  when  the  first  half  of  the  64  x  64 
Hadamard  codes  are  used  and  when  J  <  16  users. 

Many  options  exist  for  an  incorrectly  chosen  combination  of  Hadamard  codes 
shown  in  Fig.  3.2  with  Lc  =  9  and  J  =  8  users,  but  only  one  simplified  example 
will  be  given.  Since  M  =  16,  one  incorrect  option  is  to  use  Hadamard  code  one  and 
Hadamard  codes  nine  through  fifteen.  This  choice  clearly  violates  case  1  because  it 
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Figure  3.6:  A  channel  impulse  response  (CIR)  is  decomposed 
into  a  head  part,  peak  coefficient,  and  a  tail  part. 


causes  K  =  7.  Fig.  3.5  shows  v  constructed  when  using  Hadamard  code  one  and 
Hadamard  codes  nine  through  fifteen.  Note  that  with  J  users,  there  are  at  most 
M  —  J  null  tones,  and  fewer  if  the  codewords  are  chosen  incorrectly. 

3.1.2  DMT-FEQ.  The  receiver  uses  the  CIR  in  the  DMT-FEQ  to  identify 
and  mitigate  ICI  and  ISI  by  decomposing  a  channel  matrix,  C,  into  a  cyclic  part  and 
an  error  part,  Ccyci  and  CGrr,  respectively.  C  is  defined  as 

C  =  [Ct  Cc  Ch]  =  Ccycl  +  Cerr,  (3.2) 

where  Ct,  Cc,  Ch,  Ccyci,  and  Corr  are  as  defined  in  [20-22]  based  on  the  decomposition 
of  the  CIR  into  a  head  part,  peak  coefficient,  and  a  tail  part  as  shown  in  Fig.  3.6. 

Decomposing  C  into  Ccyci  and  Ccrr  is  due  to  the  ICI  and  ISI  producing  properties 
of  C.  ICI  is  produced  when  Cc  does  not  possess  a  cyclic  Toeplitz  structure,  ISI  is 
produced  from  the  preceeding  symbol  when  Ct  ^  0,  and  ISI  is  produced  from  the 
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following  symbol  when  7^  0.  Therefore,  Cerr  is  the  cause  of  ICI  and  ISI  and  is 
mitigated  with  the  DMT-FEQ. 

In  a  communications  system  with  M  total  carrier  frequencies,  N  used  carrier 
frequencies,  and  K  unused  carrier  frequencies,  the  DMT-FEQ  uses  selection  matrices, 
So  and  Si,  to  eliminate  nonrelevant  rows  and  columns  from  the  equalizer,  E.  Si  is 
defined  as 


Si  —  diag(so,  Si, . . . ,  %_i), 


(3.3) 


with 


and 


1,  if  carrier  is  used 
0,  if  carrier  is  not  used 


S0  =  Im-Si.  (3.4) 

The  selection  matrices  are  then  reduced  to  only  non-zero  columns,  where  So  is 
reduced  to  a  M  x  K  matrix  S0jred  and  Si  is  reduced  to  a  M  x  N  matrix  Si  red.  D  is 
the  diagonal  matrix  D  =  diag(o?o,<A,  ■  ■  ■ ,  rfjw-i),  where  di  =  C  (e-^"  j.  Proceeding, 
Di,red  =  S^red  •  D  ■  Sijred.  Also,  W M  and  denote  the  DFT  and  IDFT  matrix, 
respectively,  of  size  M.  Another  selection  matrix,  Zc  red,  removes  rows  in  Cerr  that 
do  not  contain  relevant  data  and  is  used  to  define  W0ired  =  S^red  •  W m  ■  Z^red  and 
Wl,red  =  S?;red  •  WM  •  Z^red. 

The  reduction  of  E  into  Eo,red  and  Eiired  is  shown  in  Fig.  3.7  and  is  performed 
in  order  to  enhance  the  speed  of  the  DMT-FEQ  by  eliminating  non-relevant  compu¬ 
tations  involving  null-tones.  Since  equalization  of  a  null-tone  will  effectively  result  in 
a  null-tone,  the  DMT-FEQ  disregards  the  null-tone  information  by  using  the  reduced 


25 


c 

(/) 

CD 

Q. 

|  used 

unused 

used 

unused 

used 

unused 

used 

ira 

unused 

.... 

used 

unused 

i 

i 

used 

: 

□ 

... 

■ 

i  i 
i  i 

unused 

j 

used 

i  n  r 

■ 

E 


l  I  l  I  l  I 

:mniij 

i  i  i  i  i  i 

TTitft" 

Eo 


-0,  red 


Figure  3.7:  An  equalizer  is  reduced  to  only  relevant  data. 

Eo,red  contains  null-tone  data  and  Ei,red  contains  used  carrier 
data. 


forms  of  the  equalizer,  E0,red  and  Ei,red-  As  described  in  [20-22]  and  illustrated  in 
Fig.  3.7,  E0,red  and  Ei,red  are  defined  as 


E, 


0  ,red 


E l^red  *  l,red  * 


t 

0  ,red’ 


(3.5) 


and 


El, red  Si, red  ‘  ^l%red  '  ®l,red> 


(3.6) 


where  the  pseudoinverse  Wjred  of  W0,red  is  defined  as 


wU=(w;i^w0i,ol)-1wl 


0  .red 


(3.7) 


and  WJ  d  denotes  the  complex  conjugate  of  W0, red- 
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Figure  3.8:  The  DMT-FEQ  removes  the  K  null-tones  from  the 
received  signal.  This  causes  jF{y}  to  be  a  vector  with  dimensions 
Nx  1,  which  will  cause  y  to  be  an  incorrect  representation  of  the 
transmission  vector,  v.  The  system  shown  above  has  null-tones 
occupying  even  numbered  carrier  frequencies. 

3.1.3  Post-Equalization  Null-Tone  Insertion.  Due  to  the  DMT-FEQ  prop¬ 
erties  shown  in  Fig.  3.7  [20-22],  a  vector  reduction  occurs  in  the  post-equalization 
JF{ y},  where  the  K  null-tones  once  present  in  the  frequency-domain  transmission 
vector,  JF{v},  are  removed  from  E{y}.  The  impact  of  the  post-equalization  null-tone 
removal  in  jF{y}  will  produce  an  inaccurate  time-domain  received  signal,  y.  This  is 
illustrated  in  Fig.  3.8. 

One  solution  that  will  correct  the  discrepancy  between  v  and  the  post  DMT- 
FEQ  y  is  to  insert  null-tones  into  jF{y}.  Of  course,  the  null-tone  insertion  needs  to 
correspond  to  the  null-tones  that  are  present  in  JF{v}.  The  null-tone  insertion,  or 
post-equalization  expansion  of  the  received  signal,  will  increase  1F{ y}  from  an  IV  x  1 
vector  to  a  M  x  1  vector.  This  is  shown  in  Fig.  3.9. 


3.1.4  Scalar  Maximum- Likelihood  Estimation.  The  maximum  likelihood  es¬ 
timate  (MLE)  of  a  received  symbol  as  described  in  [11]  was  used  in  order  to  determine 
the  optimal  estimate  of  a  received  symbol  corrupted  by  noise. 
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Figure  3.9:  Inserting  null-tones  in  jF{y}  corresponding  to 

null-tones  that  occupy  JF{v}  creates  an  equalized  signal,  y,  that 
more  accurately  represents  the  transmission  vector,  v. 


Because  the  DMT-FEQ  described  in  [20-22]  completely  eliminates  ISI  and  ICI 
due  to  transmission  through  a  non-ideal  channel,  the  post-DMT-FEQ  received  symbol 
vector,  y,  at  the  receiver  becomes 


y  =  v  +  n, 


(3.8) 


where  v  is  as  defined  in  (3.1)  and  n  is  the  processed  additive  noise  (not  necessarily 
white)  plus  residual  ISI  vector. 


In  a  communication  system  with  M  carrier  frequencies,  n  can  be  modeled  as 


n~  M{  0  ,  K  ), 


(3.9) 


where 


0 


0i 

02 


0  M 


and 


=  E[ n  nH]  ~  nn- 

n=  1 

The  noise  covariance  matrix,  K,  can  be  obtained  through  a  training  period  by 
solving  for  n  in  (3.8). 

The  distribution  of  y  then  becomes 


all 

a12  ■ 

2 

•  •  alM 

^21 

■ 

2 

•  •  a2M 

2 

aM  1 

2 

aM  2 

2 

•  •  aMM 

y  ~  A/- 1  ^  unhn  ,  K 

\  n= 1 

The  probability  distribution  function  (PDF)  of  y  given  v  then  becomes 


(3.10) 


p(y|v)  = 


(2tt)j|  Kj 


r(y-En=i“»h»)  K  1  (y-E»=i%h») 


(3.11) 


The  MLE  of  the  transmitted  symbol  of  interest  as  described  in  [11]  is  defined 


as 


uk,ML( y|v)  =  arg  max  p(y|v).  (3.12) 

uk 

Since  maximizing  p(y|v)  is  equivalent  to  maximizing  In  p(y|v),  equation  (3.12)  can 
be  rewritten  as 


Mfc,Mi(y|v)  =  arg  max  lnp(y|v).  (3.13) 

Uk 

(3.13)  reduces  to 
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Mfc,Mi(y|v) 


=  arg  max  In 

uk  [Zn)J  |  Jx| 


j 


(y  -  «nhn)H  K  1  (y  -  ^  unhn) , 


72=1 


72—1 


H 


=  arg  max  -(y  -  ^  unhn)H  K  1  (y  -  ^  ttnhn), 

uk 

72—  1  72—1 

=  arg  max  -y^K_1y  +  yflK_1  it„hn  j  +  ^ y  unhn  )  K-1y 

^  ^  ^72^72  J  K  I  ^  ^  'Unhn  J  , 


s.  72— 1 


^72—1 


s.  72— 1 


=  arg  max  -y^K  xy  +  yHK  1  ukhk  +  V\  nhr) 

uk  \ 


n^k 


ukhk  +  ^  unhn  J  K  V-  ukhk  +  ^  unhn  J  K  1  |  ukhk  +  ^  unhn  J  , 
n^k  )  \  n^k  J  \  n^k  / 

arg  max  -yHK^y  +  yIiK~1ukhk  +  yHK^  |^tinhn)  +  (ukhk)H  K-1y 


uk 

J 


\n^k 


H 


y  «nh„  J  K  V  -  (ttfchfc)"  K  1  (ufchfe)  -  K  1  I  ^  unhr, 

/  \n^k 

J 


H 


J 


H 


J 


^  J  K  uklcik  |  ^  ^  unlcin  |  K  f  ^  '  uniin 

\n^k  /  \n^k  /  \n^k 


(3.14) 


A  necessary  condition  for  the  MLE  is  to  set  the  derivative  of  (3.14)  equal  to 
zero,  which  becomes 


_d_ 

dui 


In  p(y|v) 


«fc=“fc,Mi(y|v) 


=  0. 


(3.15) 


In  order  to  clean  up  the  derivation  of  uk,ML  at  the  receiver  we  will  proceed  as 


follows:  uk  is  a  scalar,  therefore  it  is  a  term  that  can  be  factored  out,  and,  since 


duk 
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*  Jg 


of  non-itfc  terms  will  inherently  force  them  to  zero,  we  will  ignore  them.  Also,  is 
real,  so  h£  =  h*,,  where  hjjl  denotes  the  complex  conjugate  of  h^.  Proceeding  with 
these  substitutions  results  in  (3.15)  becoming 


_d_ 

duu 


yHK  1ukhk  +  (ukhk)H  K  1y  —  (ukhk)n  K  1  (ukhk) 

-  (uk hk)H  K  1  I  ^  unhn  j  -  I  unhn  |  K~lukhk  j  =  0 

\n^k  /  \n^k  /  ' 


H  TVr-1 


_d_ 

duu 


U k 


(y^K  [hk)  +  u*k  (hlK  xy)  -  u*kuk  (h^K  ^ 


uh 


hpV1  (E“”h" 

\n^k 


^ k  |  |  ^  ^  J  K  h/e 

\n^k 


=  0  (3.16) 


In  order  to  reduce  the  clutter  of  (3.16),  we  will  define  ak  =  y^K  1hfc,  a*k  = 
hfcK_1y,  bk  =  hpC_1hfc,  ck  =  h^KA1  fe^Unbn),  and 


H 


=  \J2n^kun^n)  K  1hfc.  (3.16)  becomes 


d 

dui 


{ukCLk  ~t~  ^k^k  'U'k'U'kbk  ^k^k  ^fc^fc)  0 


(3.17) 


Since  uk  =  ( ukRe  +  j  uklrn),  (3.16)  becomes  a  derivative  operation  with  respect 
to  a  complex  parameter  [4],  which  by  definition  is 


d  In  p(y|v) 
duk 


d  In  p(y|v) 


du 


_  .  d  In  p(y|v)  j  =  Q 


kRe 


du , 


kj-n 


(3.18) 


Solving  (3.18) 
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2  CJU^  ^  (UkRe  +  3ukIrn)  ak  +  ( ukRe  jukIrn )  CLk  {UkRe  jukIrn)  ( ukRe  +  3ukIm)  bk 

~  ( UkRe  —  JUklm)  Cfc  —  (UkRe  +  3ukIrn )  —  J  2  ^  +  3ukIm )  ak 

+  (UkRe  ~  3ukIm)  ak  ~  (UkRe  ~  3ukIrn )  (Mfcfle  +  3ukIrn )  4  ~  (UkRe  ~  3ukIm)  Ck 

~  (ukRe  +  jtifc,  J  4)  =  0 

2  (Tifc  i^kn^k  +  jUkImCLk  +  UkReCLk  —  jUkImak  —  UkRbk  —  jUkReUkIrnbk  +  jukReUkIrnbk 


d  UkImbk  ukReck  +  3ukIrnck  ukRf,c, 


' HeCk  3ukIrnck 

.  1  d  ( 

J  g  (  ' Ufcfle°fc  + 


■Reak  3Ukimak  UkR<pk  jukReukImbk 


Imbk  +  jcfc  -  j'4)  =  0 


+  jukReukIrnbk  +  .7  UkImbk  ukReck  +  jukImck  UkReck  3ukImckj  0 

2  (flfc  +  at  -  2 “fcflA  -  Cfc  -  4)  -  (jak  -  ja*k  -  2uklmbk  +  jt 

K  +  <)  (cfc  +  4)  t  i2afc  ,  i24  ,  ^  J2cfc  ,  i24  _  n 

- 2 - 2 - UkRebk  2 1  2 ^  ^ Ukl™k  2 — 1  “  —  U 

/  n  ,  (afe  —  at)  (cfe  -  4) 

(UkRe  —  JUkIm )  bk  —  ttkRe  —  CkRe  H  “  1  o 


_  3  uk im )  bk  —  CLkRe  -  CkRe  H  -  I  7 

(afcfle  —  CkRe  +  3akIrn  +  3ckIrn ) 

UkRe  -  3ukIm  =  - 7 - 

Ok 


bk 

~  _  (®fc  O’k) pie  J1  (®fc  “1“  Cfc) 

'Ok, ML  =  - 7 - 

bk 

K  -  Cfc) 

Ok, ML  ~  - 7 - 


7m 


h^K-1 


Ok, ML  ~ 


y  ^nhr 


(3.19) 


(3.19)  states  that  Uk,ML  is  the  optimal  scalar  MLE  of  an  individual  user’s  mod¬ 
ulated  data.  Note  that  ck  =  h^KN1  {^)Jn^kurMn^ ,  which  requires  a-priori  knowledge 
of  all  the  other  users’  4-QAM  data  and  will  be  addressed  in  section  3. 1.4. 4. 
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Output 


_  J{ynz} 

®-n-  =  FL-J{ynz}+FR0=FL-J{ynz} 

J{yz} 


=  nonzero  data 


Figure  3.10:  The  output  of  the  DMT-FEQ  inherently  con¬ 

tains  null-tones,  denoted  as  IF{y z},  which  are  passed  on  to  the 
received  signal  vector,  y.  The  null-tones  are  shown  grouped  for 
simplicity. 


3. 1.4-1  Non- Invertible  Noise  Covariance  Matrix.  (3.19)  requires  the 
noise  covariance  matrix,  K,  to  be  invertible,  which  is  the  case  if  K  is  of  full  rank, 
as  described  in  [17].  This  section  describes  the  effects  of  null-tones  on  K  and  leads 
into  section  3. 1.4. 2,  which  provides  an  alternative  approach  to  obtain  an  invertible  K 
when  presented  with  a  non-invertible  K,  and  will  lead  to  obtaining  Uk,ML- 

Since  the  equalized  signal,  y,  and  the  transmission  vector,  v,  are  known  during 
a  training  period,  the  noise  vector,  n,  is  obtained  by  solving  (3.8).  Fig.  3.10  is  a 
simplihed  representation  of  null-tones  in  the  output  of  the  DMT-FEQ,  where 


y  =  FL  •  IF{ynz\.  (3.20) 

In  (3.20),  Fl  denotes  the  left  portion  of  the  IDFT  matrix  and  IF{ ynz}  denotes  the 
frequency-domain  non-zero  data  of  the  output  of  the  DMT-FEQ.  Solving  (3.8)  will 
lead  to  the  fact  that  null-tones  that  are  present  in  y  will  also  be  present  in  n.  Setting 
the  signal  component  of  y  to  zero  and  solving  for  n  in  (3.20)  results  in 

n  =  FL  •  T{nnz\.  (3.21) 

Once  n  is  obtained,  K  becomes 
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0  =  null-tone 
□  =  non-zero  data 

Figure  3.11:  Null-tones  affect  the  received  signal  vector  y  by 
causing  linear  dependence  of  its  terms. 


K  =  E[n-nH }  =  E[FLE{nnz}  •  E{nnz}HF*}  =  FL  •  K  •  Ff ,  (3.22) 

where  K  is  an  N  x  N  noise  covariance  matrix  consisting  of  only  noise  data  from  the 
used  carrier  frequencies  and  N  denotes  the  number  of  used  carrier  frequencies.  The 
dimensions  of  are  M  x  N,  where  M  denotes  the  total  number  of  possible  carrier 
frequencies.  Therefore,  according  to  (3.22),  rank  (K)  <  N  because  rank  =  N. 
Since  N  =  M  —  K,  K  does  not  have  full  rank  which  makes  K  non-invertible. 

3. 1-4-2  Invertible,  Reduced  Noise  Covariance  Matrix.  Fig.  3.11  shows 
a  system  with  N  =  K,  each  null-tone  occupying  the  odd  carrier  frequencies,  and 
M  =  2N.  This  is  the  case  when  the  last  half  of  an  8  x  8  Hadamard  matrix  is  used  to 
construct  the  transmission  vector,  v,  when  the  users,  J  <  4.  When  the  first  half  of 
an  8  x  8  Hadamard  matrix  is  used  when  the  users,  J  =  4,  the  mill-tones  will  occupy 
the  even  carrier  frequencies.  Please  refer  to  the  frequency-domain  characteristics  of  a 
Hadamard  matrix  shown  in  Fig.  3.2  and  Fig.  3.3. 

Fig.  3.11  demonstrates  null-tones  in  JP{y}  create  linearly  dependent  terms  in 
y,  where  y\-N  =  —  y(jv+i):M-  Due  to  (3.8),  this  also  makes  ni:jv  =  —  npv+i):M-  Since 
K  =  E[ n  nT],  a  system  similar  to  Fig.  3.11  with  N  =  K ,  n  will  force  K  to  be 


34 


Figure  3.12:  As  shown  in  (3.20)  and  (3.21),  null-tones  di¬ 

minish  the  rank  of  the  noise  covariance  matrix,  K,  making  it 
non-invertible.  K  shown  is  a  result  of  null-tones  depicted  in 
Fig.  3.11,  where  N  =  K  =  A. 


composed  of  four  equal  sub-matrices,  as  shown  in  Fig.  3.12,  which  is  non-invertible 
due  to  its  rank  [17].  One  option  to  make  an  invertible  K  out  of  a  non-invertible  K 
is  to  use  a  reduced  form  of  K,  that  is  K red.  K red  is  also  defined  by  (3.22),  but  will 
originate  with  n  using  a  reduced  form,  nred.  nred  is  defined  as 


nred  =  ni  :N-  (3.23) 

Using  nred  produces  K red  that  is  defined  as 

Kred  =  E[nred  n?ed],  (3.24) 

where  K.red  is  a  JV  x  iV  matrix  with  a  rank  =  N.  Therefore,  K red  is  invertible. 

3.1. 4-3  Reduced  Received  Signal  Vector  and  Reduced  Hadamard  Codes. 
The  utilization  of  an  invertible  K.red  in  (3.19)  forces  the  requirement  of  reduced  forms 
for  the  received  signal  vector,  y,  and  individual  Hadamard  codes,  h„ .  Fig.  3.11  depicts 
an  equalized,  time-domain  received  signal,  y,  multiplied  with  a  Hadamard  code,  hn. 
Since  the  multiplication  of  the  first  halves  of  both  y  and  hn  produce  the  same  result 
as  the  multiplication  of  the  last  halves  of  both  y  and  hn,  the  reduced  terms  yred  and 
h ntred  are  respectfully  defined  as 
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Yred  yi:N, 


(3.25) 


and 


h  n,red  (3.26) 

Now  we  shall  define  a*k  red  =  h^redK“^yred,  bk)Ted  =  hlredK~^dhk>red,  and  ckn,red  = 
hjfc,redK~^  UjJOn^edj  •  Substituting  the  ak,  bk,  and  ckn  terms  from  (3.19)  with 

ak,redi  bk,red,  and  ckn,red,  respectively,  results  in 


Uk,ML 


(®Jk,red 


Ckn, red) 


^k,red^redy red  ^k^ed^red  (En/fcW^ 


n.red 


' k,red 


Vi  T 

1Afc,redxvred11^  -fed 


hT  K1 

_  ^k^ed^red 

y  red  ^2n^kUn^n  ,red  ^ 

k’ML  lt  'R'”1n 

^kjred^red^kjred 

(3.27) 


(3.27)  is  the  optimal  estimate  of  an  individual  user’s  transmitted  data,  with 
respect  to  MLE  techniques.  Note  that  ukjML  requires  a-priori  knowledge  of  all  the 
other  users’  4-QAM  data,  un^k,  which  will  be  addressed  in  section  3. 1.4. 4. 

3. 1-4-4  Scalar  Approximate  Maximum  Likelihood  Estimation.  The 
other  users’  4-QAM  data,  un^k,  in  (3.27)  requires  a-priori  knowledge  of  -un)n^fc, 
which  is  not  possible  to  estimate  in  the  maximum-likelihood  (ML)  sense  because  a 
ML  estimation  (MLE)  of  un^k  will  also  require  a-priori  knowledge  of  the  other  users’ 
4-QAM  data.  Therefore,  an  initial  approximate  MLE  must  be  acquired  for  un^k, 
which  will  be  denoted  as  un  and  will  be  defined  as 


U n 


UT  T — 1 
lln,red*^~redJ  rec ^ 

i,T  ry -1 1, 

1Ln  ,redrs~redll'nXed 


(3.28) 
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Note  that  un  in  (3.28)  is  Uk,ML  in  (3.27),  with  the  exception  of  the  Y^n^k  «nh n,red 
term.  (3.28)  provides  an  estimate  as  close  as  possible  to  a  MLE  of  wnjn^fc  due  to  the 
lack  of  a-priori  knowledge  (essentially,  ignore  other  users,  which  are  nearly  orthogonal 
anyway).  Using  the  decision  function  defined  by  (3.39)  in  section  3.1.6,  w(un)  =  un, 
and  inserting  its  result  into  (3.27)  results  in  Uk,ML  becoming 


- -l 


hT  K 

redirect 


Uk,ML  ~ 


If  red  finhn,red 


uT  TV-1  i, 

llk  ,redr^'redllk  ,red 


(3.29) 


3.1.5  Vector  Maximum- Likelihood  Estimation.  An  alternative  estimation 
approach  is  to  estimate  all  of  the  received  data  at  once,  instead  of  estimating  the 
received  data  on  a  per  user  basis  as  described  in  (3.29).  Proceeding  with  this  approach, 
the  distribution  of  the  received  data  remains  as  defined  in  (3.10).  Therefore,  (3.11) 
still  holds  true  and  the  MLE  of  the  transmitted  vector  as  a  whole  will  be  calculated, 
which  is  defined  as 


uml(y|v)  = 


U 1 

u2 

Uj 


=  arg  max  lnp(y|v), 


which  becomes 


n=l 


UAfi(y|v)  =  arg  min  (y  -  ^  unhn)H  K  1  (y  -  ^  unhr 

n=  1 

In  (3.31),  Yhn=iun hn  can  alternatively  be  represented  as 


(3.30) 


(3.31) 
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(3.31)  becomes 


v  =  Hu  = 


hi  I  h2 


V 

H  MxJ 


U  l 
U2 

Uj 


ujxl 


(3.32) 


UAfi(y|v) 


arg  min  (y  —  Hu)  K  (y  —  Hu) 

u  - v - - - ' 

L 


arg  min  L. 

U 


(3.33) 


In  order  to  determine  the  minimum  of  (3.33)  we  will  set  the  gradient  equal  to 

zero 


VUL  =  0.  (3.34) 

Since  H  is  real,  H*  =  H,  and  (3.34)  becomes 

Vu  (y hK-V  -  y^K^Hu  -  (Hu^K^y  -  (Hu^K  'Hu)  =  0, 

Vu  (yHK_1y  -  yffK_1Hu  -  uHHTK_1y  -  uhHtK_1Hu)  =  0.  (3.35) 

(3.35)  can  be  cleaned  up  by  defining  d 11  =  y^K_1H,  d  =  HTK_1y,  and  e  = 
HTK_1H.  Proceeding  with  these  definitions  [1,4] 
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(3.36) 


Vu  (— d^u  —  u^d  —  utfeu)  =  0, 
0  —  2d  —  2eu  =  0, 

Uml  =  e_1d, 


U  ML  =  (HtK_1H)  1  (HTK-1y) 


uMl  in  (3.36)  is  the  vector  MLE  of  all  users’  modulated  data  and  requires  K 
to  be  invertible,  therefore  we  will  proceed  with  reduced  forms  of  H.  K,  and  y  as  in 
section  3. 1.4. 2  and  section  3. 1.4.3,  resulting  in 


(3.37) 


3.1.6  Post-Estimation  Decision  Function.  The  output  of  the  receiver’s  es¬ 
timators,  un,  Uk,Mh  and  u Ml  will  not  be  exact  4-QAM  symbols  (±1  ±  j)  because  of 
signal  degradation  due  to  noise.  Therefore,  directly  proceeding  each  of  the  system’s 
three  estimators,  un,  Uk,ML ,  and  iik,ML,  Bayesian  hypothesis  tests  [11]  will  determine 
exact  4-QAM  symbols  based  on  the  test  statistics,  z  or  z,  where  z  is  either  un  or 
Uk,ML ,  and  the  vector  estimation  test  statistic,  z,  is  u k,ML- 

Since  the  four  values  of  the  possible  4-QAM  symbols  transmitted  are  ±1  ±  j, 
the  hypothesis  testing  for  the  test  statistic  will  be  either  +1  or  -1  for  both,  real  and 
imaginary  parts.  Therefore,  the  likelihood-ratio  test  becomes 


Ho 

z  ^  0. 


(3.38) 


H 1 
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Figure  3.13:  The  angle  between  the  real  and  imaginary  parts 
of  z  determine  9,  and  6  is  used  to  determine  the  received  4-QAM 
symbol. 

(3.38)  uses  the  real  and  imaginary  axes  as  hypotheses  decision  boundaries  to 
determine  the  Bayesian  optimal  4-QAM  symbol  received.  The  real  and  imaginary 
axes  as  hypotheses  decision  boundaries  can  also  be  interpreted  as  the  angle  between 
the  real  and  imaginary  parts  of  z.  Therefore,  a  demodulator  similar  to  [14]  was  used 
to  determine  the  received  4-QAM  symbol,  i.e.,  sign  (Real  {A}}  +  j  sign{Imag{A}}. 

Since  un,  Uk,ML ,  and  u k,ML  are  complex,  the  angle  between  the  real  and  imag¬ 
inary  parts  will  be  defined  as  9,  and  will  be  used  to  decide  if  un,  Uk,MLi  and  each 
component  of  u k,ML  are  ±1  ±  j.  As  shown  in  Fig.  3.13,  the  output  of  the  4-QAM 
decision  function,  w,  when  z  =  un  is  entered  produces 

w  (un)  =  un.  (3.39) 

Also,  the  output  of  the  4-QAM  decision  function  when  z  =  Uk,ML  is  entered  produces 

W  (■ Uk,ML )  =  Uk,ML •  (3.40) 
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Figure  3.14:  A  complete  view  of  the  TDCS  model. 


Finally,  the  decision  function  produces  the  following  result  when  (3.36)  is  entered 
to  obtain 


w  (u ml)  —  Uml, 


(3.41) 


where  the  decision  function  w( z)  in  (3.41)  acts  on  the  vector  vlml  element  by  element. 


3.2  TDCS  with  DMT-FEQ 

Fig.  3.14  illustrates  a  TDCS  system  model  that  has  inherent  null-tones  which 
can  be  exploited  by  using  DMT-FEQ.  Several  TDCS  models  were  used  in  conjunction 
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Spectrum  to  Avoid  Spectrum  to  Avoid  Forced  Null-Tones 


(a)  (b) 

Figure  3.15:  Two  different  types  of  TDCS  spectral  masks,  (a) 

A  traditional  TDCS  spectral  mask  containing  consecutive  null- 
tones  only  in  the  frequency  spectrum  to  avoid,  (b)  The  modified 
spectral  mask  contains  consecutive  null-tones  in  the  frequency 
spectrum  to  avoid  and  additional  forced  null-tones  to  maintain 
compatibility  with  DMT-FEQ  requirements. 

with  DMT-FEQ  to  try  and  identify  causes  of  incompatibility  between  DMT-FEQ  and 
a  traditional  TDCS.  The  TDCS  generally  possesses  all  pseudo-random  phases  and  a 
spectral  mask  that  only  contains  null-tones  in  the  frequency  spectrum  to  be  avoided. 

As  shown  in  Fig.  3.14,  there  are  two  major  variables  in  a  TDCS  system  model 
that  can  be  modified  for  compatibility  purposes  to  make  a  TDCS  system  mesh  well 
with  DMT-FEQ:  phase  and  spectrum  magnitude.  Therefore,  the  phase  portion  of 
TDCS  was  modeled  two  different  ways:  conjugate-symmetric  phase  and  pseudo¬ 
random  phase.  The  spectrum  magnitude  portion  of  TDCS  was  modeled  with  a  slight 
modification  to  a  traditional  TDCS  spectral  mask  that  included  forced  null-tones.  All 
of  these  variables  are  described  below. 

3.2.1  Conjugate- Symmetric  Phase.  As  described  in  [10],  a  conjugate- symmetric 

sequence  was  produced  for  the  TDCS  phase  as  follows 
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T 


(3.42) 


9 


0  9X  02 


9m 

2 


0  9 


9*2  9* 


where  M  =  64  is  the  total  number  of  carrier  frequencies  and  9  are  uniformly  dis¬ 
tributed  PNs  ranging  over  0  to  27t.  Note  that  •—  1  of  the  possible  M  phases  are 
PN,  rather  than  the  required  M  PN  phases  noted  in  [2]. 


3.2.2  Pseudo-Random  Phase.  The  pseudo-random  TDCS  phase  was  just 
that,  random  phase  values  ranging  from  0  to  27r  for  all  the  possible  M  phases. 


3.2.3  Forced  Null-Tones  Spectrum  Mask.  The  spectral  mask  in  a  traditional 
TDCS  is  created  for  the  sole  purpose  of  avoiding  a  specific  frequency  range,  repre¬ 
sented  in  Fig.  3.15  as  the  consecutive  null-tones.  The  forced  null-tones  spectral  mask 
is  a  slight  modification  to  the  traditional  TDCS  spectral  mask  and  includes  additional 
forced  null-tones  throughout  the  frequency  spectrum  to  comply  with  the  equidistant 
null-tone  rule.  The  additional  null-tones  are  added  without  violating  the  null-tone 
rules  noted  in  Section  3.1.1  [20-22]:  Case  one  is  K  >  Lc  —  1,  Lc  —  1  <  21V, and  Case 
two  is  K  >  2 IV,  Lc  —  1  >  2 N.  Both  spectral  masks  are  shown  in  Fig.  3.15. 
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IV.  Simulation  Results 


This  chapter  describes  simulation  results  of  equalizing  the  CDMA  and  TDCS 
communications  systems  described  in  Sections  3.1  and  3.2  with  the  discrete- 
multitone  (DMT)  frequency- domain  equalizer  (FEQ)  described  in  [20-22],  To  more 
accurately  describe  the  requirement  for  channel  equalization,  this  chapter  is  parti¬ 
tioned  into  three  sections.  The  first  section  simulates  the  Hadamard  encoded  CDMA 
communications  system  transmitting/receiving  through  an  ideal  channel  without  DMT- 
FEQ,  the  second  section  simulates  the  Hadamard  encoded  CDMA  communications 
system  transmitting/receiving  through  a  non-ideal  channel  without  DMT-FEQ,  and 
the  third  section  simulates  both  communications  systems  transmitting/receiving  through 
their  respective  non-ideal  channels  with  DMT-FEQ  employed.  Because  the  require¬ 
ment  for  equalization  was  demonstrated  with  the  Hadamard  encoded  CDMA  simu¬ 
lations,  the  TDCS  simulations  were  only  conducted  through  a  non-ideal  channel  in 
conjunction  with  the  DMT-FEQ.  A  diagram  of  the  simulations  performed  in  this 
chapter  is  shown  in  Fig.  4.1. 

The  non-ideal  channel  impulse  response  (CIR)  data  was  obtained  from  finite 
impulse  response  models  of  digital  microwave  radio  CIRs  in  [5].  The  initial  length 
of  this  non-ideal  CIR  data  was  Lc  =  300.  The  channel  windowing  process  placed 
the  peak  channel  coefficient  in  the  center.  Then  16  channel  coefficients  were  retained 
before  and  after  the  peak  coefficient  for  the  Hadamard  encoded  CDMA  system  and 
three  were  retained  before  and  after  the  peak  coefficient  for  the  TDCS.  All  other 
channel  coefficients  were  then  removed.  The  shortened  non-ideal  channel  length  was 
set  to  the  threshold  of  the  null-tone  rules,  where  Lc  =  33  for  Hadamard  encoded 
CDMA  simulations  and  Lc  =  7  for  the  TDCS  simulations. 

The  Hadamard  encoded  CDMA  simulations  followed  the  system  models  de¬ 
scribed  in  Section  3.1  and  were  comprised  of  combinations  of  correctly  chosen  Hadamard 
codes  that  did  not  violate  the  null-tone  conditions  noted  in  Section  3.1.1.  The  size 
of  the  Hadamard  matrix  used  was  64  x  64.  The  noise  used  in  simulations  involving 
the  first  half  of  the  Hadamard  spreading  codes,  hi:32,  was  assumed  to  be  zero-mean 
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Figure  4.1:  A  diagram  of  simulations  conducted. 


with  unit  variance  in  order  to  construct  an  invertible  noise  covariance  matrix.  Also, 
only  h33:64  were  used  in  conjunction  with  the  vector  estimator  because  hi:32  led  to 
poor  performance  when  used  in  a  Hadamard  encoded  CDMA  communications  system 
utilizing  DMT-FEQ. 

TDCS  simulations  followed  the  system  models  described  in  Section  3.2,  where 
only  one  TDCS  user  was  using  the  system’s  available  bandwidth.  The  TDCS  sim¬ 
ulations  were  split  into  two  categories:  simulations  without  a  primary  interference 
signal  and  simulations  with  a  primary  interference  signal.  Although  a  portion  of  the 
frequency  spectrum  was  restricted  from  use,  there  were  no  primary  users  transmitting 
on  the  simulated  interference  frequencies  represented  in  Fig.  3.15(a)  as  the  consecu¬ 
tive  zeros  during  the  initial  TDCS  simulations.  These  initial  TDCS  simulations  were 
conducted  for  the  sole  purpose  of  identifying  TDCS  properties  that  worked  well  with 
DMT-FEQ.  The  simulations  with  a  primary  interference  signal  were  used  to  validate 
the  findings  from  the  initial  simulations  in  a  more  realistic  TDCS  model.  Also,  be¬ 
cause  the  requirement  for  equalization  was  demonstrated  with  the  Hadamard  encoded 
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Figure  4.2:  Ideal  CIR  for  Hadamard  simulations.  Lc  =  33. 

CDMA  simulations,  the  TDCS  simulations  were  only  conducted  through  a  non-ideal 
channel  in  conjunction  with  DMT-FEQ. 

All  simulations  had  additive  white  Gaussian  noise  (AWGN)  injected  into  the 
receiver  input,  where  the  received  signal  to  noise  ratio  (SNRrx)  ranged  from  OdB  to 
40dB.  The  Hadamard  encoded  CDMA  simulations’  SNRrx  was  in  5dB  increments 
and  the  TDCS  simulations’  SNRr,x  was  in  ldB  increments. 

All  simulations  had  each  user  transmit  and  receive  30,400  of  their  respective 
4-QAM  symbols,  or,  on  the  binary  level,  each  user  transmitted  and  received  60,800 
of  their  respective  binary  data  bits. 

4-1  Ideal  Channel  without  DMT-FEQ 

Since  most  communication  systems  operate  over  channels  that  are  non-ideal,  the 
results  in  this  section  are  used  as  a  baseline  for  the  subsequent  simulation  results  that 
model  more  realistic  environments  for  a  Hadamard  encoded  CDMA  communications 
system. 

The  Hadamard  encoded  CDMA  system  is  shown  in  Fig.  3.1,  where  the  only 
exception  was  that  the  DFT,  DMT-FEQ,  null-tone  insertion,  and  the  IDFT  were  not 
used.  The  number  of  synchronous  users,  J  =  1,2,4,8,16,  and  32.  The  ideal  channel 
impulse  response  (CIR)  in  these  simulations  is  of  length,  Lc  =  33,  as  shown  in  Fig. 
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(a)  (b) 

Figure  4.3:  P&  vs.  SNRRx  results  for  a  Hadamard  encoded 
CDMA  communications  system  with  an  ideal  CIR,  not  using 
the  DMT-FEQ,  and  using  the  scalar  estimator,  (a)  h1:32  were 
used,  (b)  h33:64  were  used. 

4.2.  Although  Lc  =  33  in  Fig.  4.2,  zero  signal  delay  and  zero  signal  distortion  will 
occur  when  a  signal  is  applied  to  the  ideal  CIR. 

4-1.1  Scalar  MLE,  hi:32.  Although  an  ideal  CIR  was  used,  the  probability 
of  bit  error,  Pb ,  versus  SNRRx  results  in  Fig.  4.3(a)  indicate  high  bit  error  rates  when 
a  CDMA  system  accommodates  16  or  more  synchronous  users  when  h1:32  was  used. 

4-1-2  Scalar  MLE ,  h33:64.  Unlike  the  Pb  versus  SRTRf{x  results  in  Fig.  4.3(a), 
the  results  in  Fig.  4.3(b)  show  Pb  improvement  as  SNRRx  increases.  Therefore,  these 
results  suggest  case  one  or  case  two  and  the  equidistant  null-tone  rules  must  be  ad¬ 
hered  to  for  proper  equalization  by  DMT-FEQ. 

4-1-3  Vector  MLE ,  h33:64.  Results  from  using  the  vector  estimation  are 
shown  in  Fig.  4.4. 
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Figure  4.4:  Pf,  vs.  SNRrx  results  for  Hadamard  encoded 

CDMA  communications  system  with  with  an  ideal  CIR,  no 
DMT-FEQ,  and  using  the  vector  estimator.  1133,64  were  used. 


4-2  Non-Ideal  Channel  without  DMT-FEQ 

Many,  if  not  all,  fielded  communications  systems  are  not  afforded  the  luxury  of 
operating  through  ideal  channels.  Therefore,  this  section  simulated  the  Hadamard  en¬ 
coded  CDMA  communications  system  operating  through  a  non-ideal  channel  without 
utilizing  DMT-FEQ  and  is  shown  in  Fig.  3.1,  where  the  only  exception  was  that  the 
DFT,  DMT-FEQ,  null-tone  insertion,  and  the  IDFT  were  not  used.  The  CIR  used 
for  simulations  is  shown  in  Fig.  4.5(a),  where  the  number  of  users,  J  =  1,2,4,8,16, 
and  32. 
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(a)  (b) 

Figure  4.5:  Non-Ideal  CIRs.  Original  length  of  CIR  was  Lc  = 
300  and  was  obtained  from  [5].  (a)  Non-ideal  CIR  used  for 

Hadamard  simulations  shortened  to  Lc  =  33  in  order  to  barely 
satisfy  the  case  one  null-tone  rule  noted  in  section  3.1.1.  (b) 
Non-ideal  CIR  used  for  TDCS  simulations  shortened  to  Lc  =  7. 


4-2.1  Scalar  MLE,  h1:32.  The  Pb  versus  SNRRx  results  shown  in  Fig.  4.6(a) 
reflect  the  consequences  of  not  equalizing  a  signal  after  transmission  through  a  non¬ 
ideal  channel.  The  results  indicate  poorer  performance  in  terms  of  bit  error  rate  when 
no  equalization  was  used  because  Pb  increased  from  an  approximate  average  of  1  in 
50,  shown  in  Fig.  4.3(a),  to  an  approximate  average  of  1  in  5,  shown  in  Fig.  4.6(a). 

4-2.2  Scalar  MLE,  h33:64.  The  Pb  versus  SNRRx  results  shown  in  Fig.  4.6(b) 
also  reflect  the  impact  of  not  equalizing  a  signal  after  transmission  through  a  non¬ 
ideal  channel.  The  results  indicate  poorer  performance  in  terms  of  bit  error  rate 
when  no  equalization  was  used  because  Pb  increased  from  the  baseline  results  shown 
in  Fig.  4.3(b)  to  an  approximate  average  for  users  J  >  4  of  1  in  5.  The  increase  in  Pb 
due  to  not  utilizing  the  DMT-FEQ  is  shown  in  Fig.  4.6(b). 


49 


(a)  (b) 

Figure  4.6:  Pb  vs.  SNRRx  results  for  a  Hadamard  encoded 
CDMA  communications  system  with  a  non-ideal  CIR,  no  DMT- 
FEQ,  and  using  the  scalar  estimator,  (a)  hi:32  were  used,  (b) 
h33:64  were  used. 

4-2.3  Vector  MLE,  h33;64.  The  Pb  versus  SNRRx  results  shown  in  Fig.  4.7 
illustrates  the  consequences  of  not  using  an  equalizer  in  a  Hadamard  encoded  CDMA 
communications  system  operating  through  a  non-ideal  channel. 

4-3  Non- Ideal  Channel  with  DMT-FEQ 

This  section  simulated  the  two  communications  systems,  Hadamard  encoded 
CDMA  and  TDCS,  operating  through  their  respective  non-ideal  channels  and  utilizing 
DMT-FEQ.  The  simulated  Hadamard  encoded  CDMA  communications  system  is 
shown  in  Fig.  3.1,  where  the  number  of  users,  J  =  1,2, 4, 8, 16,  and  32.  The  channel 
data  used  in  this  section  for  the  Hadamard  encoded  CDMA  communications  system 
was  the  same  channel  data  that  was  used  in  section  4.2,  and  is  shown  in  Fig.  4.5(a). 
The  channel  data  used  for  the  TDCS  simulations  is  shown  in  Fig.  4.5(b). 

4-3.1  Scalar  Estimator,  h1:32.  The  Pb  versus  SNRRx  results  shown  in 
Fig.  4.8(a)  demonstrate  improvement  over  the  Pb  versus  SNRRx  results  shown  in 
Fig.  4.6(a). 
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Figure  4.7:  P&  vs.  SNRRx  results  for  a  Hadamard  encoded 
CDMA  communications  system  with  a  non-ideal  CIR,  no  DMT- 
FEQ,  and  using  the  vector  estimator.  h33:64  were  used. 

4-3.2  Scalar  Estimator,  h33:64.  The  Pb  versus  SNRRx  results  shown  in 
Fig.  4.8(b)  demonstrate  improvement  for  a  Hadamard  encoded  CDMA  communica¬ 
tions  system  for  all  of  the  number  of  system  users  compared  with  the  Pb  versus  SNRRx 
results  shown  in  Fig.  4.6(b). 


4-3.3  Vector  Estimator,  h33:64.  The  results  from  using  the  vector  estimator 
are  shown  in  Fig.  4.9,  and  are  an  improvement  over  the  scalar  estimator  shown  in 
Fig.  4.8(b). 


4-3-4  Modified  TDCS  with  No  Primary  Interference.  The  simulated  TDCS 
is  shown  in  Fig.  3.14,  where  the  only  exception  here  is  that  the  primary  interference 
transmit  antenna  was  removed.  Also,  the  number  of  TDCS  users,  J  =  1. 
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(a)  (b) 

Figure  4.8:  P&  vs.  SNRRx  results  for  a  Hadamard  encoded 
CDMA  communications  system  with  with  a  non-ideal  CIR,  us¬ 
ing  DMT-FEQ,  and  using  the  scalar  estimator,  (a)  hi:32  were 
used,  (b)  h33;64  were  used. 

The  simulation  results  shown  in  Fig.  4.10  show  that  a  traditional  TDCS  does  not 
work  well  with  DMT-FEQ  because  of  the  relatively  high  bit  error  rates  compared  with 
a  modified  TDCS  with  conjugate-symmetric  phase  and  a  forced  null-tones  spectral 
mask.  Note  that  these  modifications  to  TDCS  produced  better  results  and  retains 
the  required  noise-like  appearance. 


4-3.5  Modified  TDCS  with  Primary  Interference.  The  simulated  TDCS  with 
no  primary  interference  signal  is  shown  in  Fig.  3.14,  where  the  number  of  TDCS  users, 
J  =  1. 

The  simulation  results  shown  in  Fig.  4.11  show  that  the  system’s  bit  error  rates 
remain  consistent  regardless  of  the  signal  to  interference  ratio  (SIR). 
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Figure  4.9:  Pb  vs.  SNRrx  results  for  a  Hadamard  encoded 
CDMA  communications  system  with  with  a  non-ideal  CIR,  us¬ 
ing  the  DMT-FEQ,  and  using  the  vector  estimator.  h33:64  were 
used. 


Figure  4.10:  TDCS  simulation  results  with  a  non- ideal  CIR, 
Lc  =  7,  shown  in  Fig.  4.5(b)  and  using  DMT-FEQ.  J  =  1  user 
with  no  primary  interference  signal. 
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Bit  Error, 


Figure  4.11:  Simulation  results  with  a  modified  TDCS  hav¬ 
ing  a  conjugate- symmetric  phase  and  forced  null-tones  spec¬ 
tral  mask.  Simulations  were  conducted  with  the  non-ideal  CIR 
shown  in  Fig.  4.5(b),  using  DMT-FEQ,  J  =  1  user,  and  a  pri¬ 
mary  interference  signal.  The  spectral  mask  used  is  shown  in 
Fig.  3.15(b). 
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V.  Conclusion 


Results  from  Chapter  4  suggest  the  discrete  multitone  frequency-domain  equalizer 
(DMT-FEQ)  works  best  in  a  multi-user  Hadamard  encoded  CDMA  communi¬ 
cation  system  environment  when  the  vector  estimation  is  performed,  rather  than  the 
scalar  estimation.  Also,  the  last  half  of  the  available  spreading  code  set  produced 
better  performance  than  the  first  half  of  the  spreading  code  set. 

The  DMT-FEQ  may  also  be  suitable  for  covert  applications  by  modifying  the 
traditional  TDCS  containing  all  pseudo-random  phases  to  include  only  consecutive 
null-tones  in  its  spectral  mask.  When  using  DMT-FEQ  with  TDCS,  simulation  results 
show  that  modifications  to  TDCS’  phase  component  to  conjugate-symmetric  phases 
and  modifications  to  TDCS’  spectrum  mask  component  to  include  forced  null-tones 
provide  acceptable  results.  This  may  be  a  viable  alternative  given  that  half  of  the 
required  phases  were  assigned  as  pseudo-random  numbers,  then  were  conjugated  and 
mirrored  in  the  remaining  half  of  the  required  phases,  which  may  preserve  the  required 
noise-like  appearance  of  TDCS  signals  in  the  time-domain. 
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